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1 Introduction 

Photodissociation regions (PDRs) form on the surfaces of molecular clouds when- 
ever these clouds are struck by far-ultraviolet radiation from hot young stars. 
These regions are characterized by the transition from hot, ionized gas to cold, 
molecular gas as the far-ultraviolet field is attenuated farther from the ionization 
front. The temperature profile of the PDR varies depending on the density and 
strength of the FUV field. Derivation of this profile must take into account the 
local heating and cooling, the chemical equilibrium, and radiative coupling to 
other layers within the structure. The Orion Bar is a dense molecular structure 
at the southeast boundary of the Orion Nebula. Early theoretical models of the 
Orion Bar by Tielens and Hollcnbach [1] (for a density of 2.3 x 10 4 and UV field 
strength (Go) of 10 5 ) predict temperatures of ^1000 K at Ay=0-2, dropping to 
less than 100 K by A v =4. 

Ground state rotational lines of H2 are good temperature probes of moder- 
ately hot (200-1000 K) gas. The low A-values of these lines result in low critical 
densities while ensuring that the lines are optically thin. ISO observations of H 2 
rotational lines in PDRs reveal large quantities of warm gas that are difficult to 
explain via current models[2], but the spatial resolution of ISO does not resolve 
the temperature structure of the warm gas. We present and discuss high spatial 
resolution observations of H2 rotational line emission from the Orion Bar. 



2 Observations and Data Reduction 

We mapped the H 2 v = 0-0 S(l) and S(2) lines at 17.03 \im and 12.28 um toward 
the Orion Bar in 2002 December. We made the observations using the Texas 
Echelon Cross Echelle Spectrograph (TEXES, [3]) on the 3m NASA Infrared 
Telescope Facility (IRTF). The slit width was 2.0" for the 0-0 S(l) line and 1.4" 
for the 0-0 S(2) line. The spectral resolution of our data is 5.5 km s" 1 at 17.03 
\im and 4.7 km s _1 at 12.28 |0.m. On source integration times per position were 
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160 and 80 seconds for the 0-0 S(l) and S(2) lines respectively. We oriented the 
slit parallel to the Orion Bar ionization front, at a position angle of 45°. We 
mapped by stepping the telescope from northwest to southeast in 0.7" (for S(2)) 
and 1" (for S(l)) steps to create 40" long scans. We reduced the raw images of 
cross-dispersed spectra using the standard TEXES pipeline reduction program 
[3]. We then smoothed our data to a spatial resolution of 2". The 0,0 position 
for the maps is at R.A. = 5 h 35 m 19.7 s , Dec. = -5° 25' 28.3" (J2000.0) and the 
mapped region runs from 13" northwest to 27" southeast of this position. 

3 Results 

3.1 Morphology 

Figure 1 shows the distribution of u=0-0 S(l) and S(2) intensity along with 
distribution of intensity in the 2.12 (am v = 1-0 S(l) line [4], resampled onto the 
same grid at the same spatial resolution. The dominant feature in the maps is the 
bright horizontal (northeast-southwest) ridge centered at y~-2". This feature has 
a sharp northwestern edge, rising at most positions from a low or undetectable 
level to half of its peak intensity in <2" (0.004 pc). The total thickness of the 
bright ridge is ~10" (0.02 pc). 

There is remarkable agreement in the intensity distributions, not only be- 
tween the 0-0 S(l) and 0-0 S(2) lines but also between these lines and the v=l-0 
S(l) line. The maps agree to within the uncertainties about the location and 
width of the bright ridge, and about the presence and extent of lower-level ex- 
tended emission. To the left of its center, the bright ridge contains a clumpy "V" 
structure whose overall morphology echoes in all three lines. 

3.2 Temperatures and Column Densities 

Table 1 lists the temperatures and column densities for the areas labeled in 
Figure 1. Derived excitation temperatures range from roughly 400 to 700 K, 
with most areas at around 600 K. Over the 20" range in depth where we are 
able to measure temperature, we do not see any systematic trend in temperature 
with depth into the molecular cloud (to the southeast and farther from 9 1 G). 

At most positions the temperature derived from the 0-0 S(l)/0-0 S(2) line 
intensity ratio agrees to better than 50 K with the temperature derived from the 
1-0 S(l)/0-0 S(l) line intensity ratio despite a factor of 4 difference in F, upper 
between the 0-0 S(2) and 1-0 S(l) transitions. Looking at the entire dataset, 
we see the most variability between the two derived temperatures close to the 
ionization front. Because of the small variation in temperature over the mapped 
region of the PDR, the 0-0 S(2) line intensity map roughly traces the column 
density of warm H 2 . 

3.3 Line widths and velocities 

With the spectral resolution available with TEXES, we are able to resolve the 
ground vibrational state H 2 lines. Physical linewidths (after deconvolving the 
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Fig. 1. Maps of integrated intensities for the S(l) and S(2) pure rotational lines taken 
with TEXES along with the v — 1-0 S(l) line [4]. Labelled axis are in arcseconds. The 
scans were taken perpendicular to the ionization front at a position angle of 45° . The 
0,0 position in the maps corresponds to R.A. = 5 h 35 m 19.7 s , Dec. = -5° 25' 28.3" 
(J2000.0). The wedge of the greyscale goes from the 1 sigma noise to the maximum 
value of the maps, or 0.33 to 9.6 x 10~ 4 ergs cm" 2 s" 1 sr" 1 for 0-0 S(l), 0.22 to 8.2 x 
10~ 4 ergs cm~ 2 s _1 sr _1 for 0-0 S(2), and 0.25 to 4.4 x 10~ 4 ergs cm~ 2 s _1 sr _1 for 1-0 
S(l). The middle panel labels the areas averaged together to increase signal-to-noise 
for further analysis. 

instrument profile) are 3-6 km s _1 , in agreement with lincwidths expected for 
optically thin, thcrmalized gas at 600 K. Line widths for the two transitions agree 
to within 1-2 km s _1 , with neither line being systematically wider than the other. 
We measure Vlsr f° r the H 2 lines to be 10-11 km s _1 . Our observed line widths 
and velocities are in agreement with published values for both molecular lines 
(NH 3 [5], CS, and isotopes of CO [6]) and carbon recombination lines [7]. 
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Table 1. Physical Conditions 



Observed Intensities 
Area 0-0 S(l) 0-0 S(2) 1-0 S(l) T ex a T ox b N(H 2 ) C 
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" 4 erg cm 




K 


K 


10 20 cm- 2 


A 


7.7 


6.6 


3.6 


476 


630 


7.7 


B 


3.8 


4.3 


0.49 


591 


554 


3.1 


C 


1.8 


2.1 


1.3 


616 


661 


1.4 


D 


1.1 


1.3 


0.94 


668 


673 


0.81 


E 


2.6 


1.8 


<0.2 


413 


<525 3.1 



a Excitation temperature determined from I[0-0 S(l)]/I[0-0 S(2)] assuming optically 
thin emission with o/p=3. 

b Excitation temperature determined from I[l-0 S(l)]/I[0-0 S(l)] assuming optically 
thin emission with o/p=3. 

c Column density of warm Eh for thermalized gas at the temperature determined from 
I[0-0 S(l)]/I[0-0 S(2)]. 



4 Discussion 

PDR models predict that H 2 emission arises at a depth in the PDR with steep 
temperature gradients [1,2]. We do not see evidence for such spatial gradients 
even though we have resolved the PDR. The agreement of temperatures derived 
from the 0-0 S(l)/0-0 S(2) and 1-0 S(l)/0-0 S(l) line intensity ratio indicates 
that the emitting gas is isothermal. We re-interpret the 400-700 K temperatures 
found for galactic PDRs in the large ISO beam [8] as emerging from isothermal 
gas, and not the beam-average of a temperature gradient. Future PDR modeling 
efforts must be able to explain the ubiquity of 600 K H 2 in PDRs. 
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